Magnetic instability of a two-dimensional Anderson non-Fermi Liquid 
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We show that in the Anderson model for a two-dimensional non- Fermi liquid a magnetic instability 
can lead to the itinerant electron ferromagnetism. The critical temperature and the susceptibility of 
the paramagnetic phase have been analytically calculated. The usual Fermi behaviour is re-obtained 
taking the anomalous exponent to be zero. 

PACS numbers: 



Since the discovery of high-Tc superconductors there 
has been a great interest in the study of the non-Fermi 
liquid model in two. jdimensions (2D). The model pro- 
posed by AndersonBfl has a phenomenological charac- 
ter because the form of the Green's function has been 
taken from the one-dimensional Luttinger liquid. Later 
the model has becBi used for both the normal and super- 
conducting stateu^, and the usual BCS critical temper- 
ature was re-obtained for the case of a zero anomalous 
exponent on the Green function. 

However, the electron-hole channel and the possibil- 
ity for the occurrence of the itinerant electron ferro- 
magnetisnL-iwas not studied for this model. A recent 
calculationia of the magnetic susceptibility of a ID sys- 
tem, done with the renormalization group method and 
the Monte Carlo simulation, showed that the magnetic 
susceptibility x(T) becomes constant as T — > 0, but with 
an infinite slope, in contrast with the normal Fermi liq- 
uid. 

In the present paper we shall assume the existence of 
a 2D non- Fermi liquid described by the Green function: 



G(k,c.) 



g{a)e 
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where a is the anomalous exponent, (p = — 7ra/2 is a 
phase factor introduced in order to respect the time re- 
versal symmetry of the Green's function, £(k) is the ki- 
netic energy of the electrons and g{a) — 7rQ;/2sin (7ra/2) 
is a factor introduced to preserve the equal time anti- 
comutation relations of the electrons. The form of the 
Green's function is assumed to be given by Eq. (|l|) as 
long as —ujc < w < ll>c, being the energy cutoff spe^. 
cific for the model. The exponent a is non-|U|iiiY|ecsal3'tl 
and is considered to has a value < a < l/2llla'QLli3. 

In order to study the electron-hole instability we cal- 
culate the polarisation n(q, oj defined by 
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[iuJn - iuj„i - C(k - q)]^ " 
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with ^(k) = k'^ /2m — fJ, + nV/2, fi being the chemical 
potential and nV/2 the Hartree energy of the n electrons. 
The summation over the Matsubara frequencies can be 
substituted by an complex integral and we obtained for 
the polarisation: 



n(q, iUjn) = 2e~ 



dx 
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Following the same way as in Ref. |ll| we calculate the 
polarisation as 



n(q,icj,„) = 2iV(0)e- 
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where 



5i(/3,/i) = 5^(-l) 



52(/3,/i) = 5J(-1) 

m=0 
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n(q,u;„) = e 
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where 



(2nr 



S{k,q,iuJm) (2) 



m— 



with p, = fi — nV 12. 



Tl + 1)]2"' 



(8) 
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The paramagnetic susceptibility given by 

x(T) = -2Aigi?en(0,0) (9) 
has been calculated as 

X(T) = -2/ioiV(0)^5- — ^cos(7ra)— — —Sl{f3,^i) 



r(l - a) 



(10) 



where N(0) is the density of states. In the limit a ^ 
we calculated from Eq. (0) 



S',{f3,fl) = fF^D{nV/2-^^) 



(11) 



where fp^oix) is the Fermi-Dirac distribution function, 
and using this result the paramagnetic susceptibility at 
T ^ becomes 



(12) 



a result which is in fact the Pauli paramagnetic suscep- 
tibility. 

The transition temperature can be calculated by con- 
sidering also the equation for the particle density in order 
to eliminate the effective chemical potential fl. The gen- 
eral equation for the electron density is 



(2^)2 



G(k, iuJn) 
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where G(k, ia;„) is given by Eq. (|l|). Performing the 
summation over the Matsubara frequencies w„ and the 
integral for a constant density of states we obtain 

(rp N{0)g{a) sin(7ra/2) 
n[T, ^i) = r(a) S'o(P, N (14) 

If we consider Sq{(3, jl) well approximated by its value at 
a = we have 



and we get 



5o(AA) = ^ln[l + e'3'^] (15) 



where 



iV(0)g(a) ^, sin(7ra/2) 



In order to calculate the critical temperature bellow 
the itinerant ferromagnetism appear we will use the re- 
lation 



1 + VRe\l{q_, iuJn) = 



(17) 



in the limit q = and icun = 0. We mention that accord- 
ing to the Mermin- Wagner theorem the ferromagnetic 



phase is destroyed by the fluctuations and Tc is in fact 
zero. However, the weak three dimensionality which is 
present in High- Tc materials gives us the possibility to use 
Eq. (|r^ ) for the calculation of the critical temperature. 
We can introduce this effect phcnomenologically taking 
instead of Eq. (0) the equation 1 + VRen{0, 0) = 0.01 
which corresponds to 10^ times enhancement of param- 
agnetic susceptibility. The qualitative phase diagram is 
not affected by these details. We calculated from Eq. (Q) 
the real part of H(q, iw„) as 



i?en(0,0) = 27V(0)r(2a- 1) 



g^{a) sin7r(l — a) 



(18) 



If we consider again that Si{j3,jl) is well approximated 
by its value for a ^ we have 



p0ii 
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which together with Eq. ( |l8| ) gives for the real part of 
the polarisation the value 



i?eH(0,0) = A{a] 



e/3/i + 1 
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with 



A{a) = 2r(2a- 1) 



5'(«) 
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Using this result we calculate the critical temperature as 
n 1 



B{a) 



1 - 



V\A(a)\ 
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From Eq. (|2^) we get a critical condition for the existence 
of the ferromagnetic order expressed as 



K < 



1 



\A{a)\ 



(22) 



In the limit a — > the critical temperature becomes iden- 
tical with the one obtained for a three dimensional itin- 
erant electron ferromagnet. This result is given by the 
fact that in both cases the integral over the energy vari- 
able is performed at the Fermi surface using a constant 
density of states. For a 2D electron system a more real- 
istic description will be the one in terms of a van Hove 
density of state. We expect as for the superconducting 
critical temperature calculated in Ref. ^ that the ferro- 
magnetic critical temperature will be enhanced by the 
energy dependence of the density of states. Anyway such 
a calculation is much more difficult and the results will 
be published in another work. 
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